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ABSTRACT 

It is pointed out that in doing the factorial moment analysis with non- 
integer partition M of phase space, the influence of the phase-space varia- 
tion of two- (or more-) particle correlations has to be considered carefully. 
In this paper this problem is studied and a systematic method is devel- 
oped to minimize this influence. The efficiency and self-consistency of this 
method are shown using the data of 250 GeV/c 7r + p and K + p collisions 
from the NA22 experiment as example. 
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1 Introduction 

Since the pioneer work of Bialas and Peschanski' 1 ' , the anomalous scaling of factorial 
moments (FM) defined as 

F / A x 1 ^ {n m (n m - 1) ■ ■ ■ (n m - q + 1)) 

has been searched extensively in high energy experiments' 2 ' . In the above equation, M 
is the partition number of a phase space region A in consideration, 5 = A/M is the 
size of a sub-cell, n m is the number of (charged) particles falling into the mth sub-cell. 

After a long period of debate, it has been recognized that in studying the anomalous 
scaling of FM the anisotropy of phase space should be taken into account' 3 ', i.e. if 
dynamical fluctuation does exist in high energy hadron-hadron collisions then it should 
be anisotropic and the corresponding fractal is "self-affme" W . This means that the 
anomalous scaling of FM can be observed when and only when the phase space is 
divided in an appropriate anisotropical way. 

Let the three phase space variables be denoted by p a , pb, p c , the corresponding shrink- 
ing ratios are A a , A& and A c respectively. The anisotropy (self- affinity) of dynamical 
fluctuation can be charcterized by the so-called roughness or Hurst exponents' 4 ' 

In A 

Hij = hiX' (^'i = a ' & or a ' c or & ' c )> ( 2 ) 

with 

Xi < \j, < < 1. (3) 

These exponents can be deduced from the experimental data by fitting three one- 
dimensional second-order factorial-moments to the saturation curve' 5 '. If self-affine 
fluctuations of multiplicity do exist in multiparticle production, exact scaling , i.e. 
a straight line in ln-FM versus InM, should be observed when and only when the 
three-dimensional analysis is performed with the phase space divided anisotropically 
according to the value of these Hurst exponents. 
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In order to confront this assertion with experiment, an important point has to be 
noticed, i.e. when H ^ 1 the partition numbers M; L in the three phase space directions, 
in general, cannot be integer numbers simultaneously. 

It turns out that the case of 250 GeV/c 7r + p and K + p collisions from NA22 exper- 
iments is special. The Hurst exponents determined from the one-dimensional second- 
order factorial-moments in these collisions are' 6 ] H yPt = 0.48 ± 0.06 ; H yif = 0.47 ± 
0.06 ; H PtLp = 0.99 ± 0.01 . These Hurst exponents can be simply approximated by 
Hyj = 1/2, (j = p t , ip), and H pttp = 1. In this approximation, the following integer val- 
ues for Mi can be used simultaneously in higher-dimensional analysis: M y = 1, 2, 3, . . . 
and M Pt = M v = 1,4,9, . . .. 

In Ref. [6] this approximation has been used and the corresponding 3-D analysis 
has been preformed. The self-affine 3D lnF 2 vs. InM thus obtained has been fitted to 
a straight line and the result is basically confirmative. However, it can be seen from 
the above that the restriction of M to integer values results in a 3-D plot with only 
7 points. Omitting the first point to get rid of the influence of transverse momentum 
conservation' 7 ! only 6 points remain, which are unfortunately insufficient for a precise 
check of the linearity of the plot. 

In order to get more points in the self-affine 3-D ln-FM versus InM plot, non- 
integer values of M have to be used. In so doing an important problem arises, i.e. the 
influence of the variation in phase space of two- (or more-) particle correlations has to 
be considered carefully. In this paper we will take this problem into account and try to 
develop a systematic method in minimizing the influence of the deviation of two- (or 
more-) particle correlations from a constant behaviour. 

In section 2 the definition of factorial moments with non-integer partition number 
M is introduced and the problem of how to use it in real experimental data analysis is 
discussed. It is shown that, due to the variation of two-particle correlations, a correction 
factor is needed. How to extract this factor from experimental data is discussed in 
sections 3 and 4. Concluding remarks are given in section 5. 
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2 Non-integer FM analysis 

For definiteness, let us consider the second order FM in one-dimensional phase space, 
e.g. rapidity y: 

F 2 (6y) = F 2 (M) = iE {nm ^ m - 1) \ (4) 

1V1 m =l \'hn/ 

where M = Ay/Sy . 

In the ideal case, factorial moments F q only depend on the bin width 5y, but not on 
the position of the bin on the rapidity axis. If that is the case, the result of averaging 
over all the M bins as in Eq. (4) is equal to that of averaging over N bins with N < M. 
That is, F 2 (M) = F 2 {N, M), where 

F 2 (N,M) = ^Y, <Wm( / Wm " 1)> , (N = M-a, < a < 1). (5) 

iV m =l \ n m) 

This equation can be used as the definition of an FM for any real value of M®. 

As is well known, even in the central region the rapidity distribution is not flat. The 
shape of this distribution influences the scaling behaviour of FM. At first, Fialkowski^ 
has suggested to use a correction factor to minimize this influence. Later, the cumulant 
variable 

/ \ Ilp(y) d y f ^ 

x(v) = (6) 

was introduced' 10 ^ , which has a flat distribution and the Fialkowski method of correction 
factor has been substituted. It turns out, however, when the FM analysis with non- 
integer partition is concerned, a similar factor has to be introduced again. 

To see this, let A denote the phase space region in consideration, 5 m the mth bin, 
Pi(yi) and p 2 (yi, y 2 ) the one- and two-particle distribution functions, respectively. Then 
we have 

(n) f 

pi(yi) = t~( — ru / P 2 (yuy2)dy 2 ; (7a) 
(n(n — 1)) ja 

(n m ) = I pi(y)dy = t^-ttt / dy 2 [ dy 1 p 2 (y 1 ,y 2 ); (76) 
Js m {n(n - 1)) J a Js m 

(n m (n m - 1)) = / dy 2 dy x p 2 {y x ,y 2 ). (7c) 



After transforming to the cumulant variable, (n m ) becomes a constant, independent 
of m. Comparing the second and third of the above equations, it can be seen clearly that 
due to the difference in the integration region over y 2 , even though (n m ) is constant, 
(n m (n m — 1)) is in general not constant. 

As an example, in Fig.l are shown /^t^) = ( n m( n m — 1))/ ( n m) 2 distributions in ra- 
pidity y transformed into the cumulant variable. Data are from 7r + p and K + p collisions 
at 250 GeV/c (NA22). The four figures correspond to M — 8,16,32,64, respectively. 
It can be seen that, although the cumulant variable has been used so that the av- 
erage distribution (n m ) is flat, the distribution of the second-order factorial moment 
/ 2 (m) = (n m (n m - 1))/ (n m ) 2 is not flat. 

Note that in the definition of F q , equations (4) and (5), a horizontal average has 
been taken. When the partition number M is an integer, the horizontal average is 
over the whole region A, cf. Eq.(4). The variation of (n m (n m — l))/(n m ) 2 is thus 
smeared out, and no correction is needed. On the contrary, when M is non-integer, 
the horizontal average is performed only over part of the region, cf. Eq.(5), and the 
influence of the variation of (n m (n m — l))/(n m ) 2 becomes essential. This influence is 
exhibited clearly in Fig.2Q In these figures, lnF 2 versus InM for y, p t and ip are shown 
with both integer (full circles) and non-integer (open circles) values of M. All these 
plots show a "sawtooth" shape. The "sawteeth" come from non-integer M. They lie 
above the smooth curve of integer M in the cases of y and p t , but reach from above the 
smooth curve of integer M to below it in the case of <p. 

In the following we will introduce a correction factor to minimize the influence of 
the variation of (n m (n m — 1))/ (n m ) 2 and remove the "sawteeth" from the lni^ versus 
InM plots for real M. 

1 As in the case of integer M, the errors of the data points are correlated. This problem has been 
discussed in Ref. [6] . 
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3 Correction factor for the (n m (n m — 1)) / (n m ) 2 distri- 
bution 

From the definition of F 2 (N,M), Eq.(5), it can be seen that when M is non-integer 
(M = N + a,0<a<l), only iV bins are included in the horizontal average, i.e. only 
r = N/M of the whole region A has been taken into account. The result is therefore 
influenced by the variation of (n m (n m — 1))/ (n m ) 2 . We try to minimize this influence 
by introducing a correction factor R(r) and define 

In order to extract the correction factor R(r) = R(N/M) from the experimental 
data, let us choose two integers N' and M', satisfying 

~ — = r. 9 

M' M K ' 

One is then ready to calculate F 2 (N',M') according to Eq.(5). 
Let us call the ratio 

^7 E (n m (n m - l))/(n m ) 2 
C(N>, M') = (10) 



ip E_(n m (n m - 1))/K Vi 



M 

Y. 

m=l 

of F 2 (N', M') to F 2 (M') correction matrix. It is the ratio of the FM averaging only over 
N'/M' of the whole region to the FM averaging over the whole region. Having Eq.(9) 
in mind, it is reasonable, therefore, to take C(N'.M') as the approximate value of the 
correction factor R(r) = R(N/M). 

In the left column of Fig. 3 are shown the correction matrix C(N, M) with iV and 
M both taking integer values as function of N/M for M — 3, 4, . . . , 40. It can be seen 
from the figures that all the points for different M's lie in a narrow band, i.e. C(N, M) 
is mainly a function of r = N/M and is only weakly dependent on M. This indicates 
that the correction factors R y (r), R pt (r) and R v (r) for non-integer M can be obtained 
from the interpolation of the corresponding C(N, M)'s. 
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It turns out that in order to eliminate the "sawteeth" in the lnF 2 versus InM plot 
appropriately, the interpolation curve should be chosen in the following way^ 11 ': When 
the "sawteeth" in the lnF 2 versus InM plot lie above the smooth curve of integer M, 
as in the case of y and p t , the upper boundary of the C(N, M) band has to be used for 
the interpolation; while when the "sawteeth" in the lnF 2 versus InM plot reach from 
above the smooth curve of integer M to below it, as in the case of </?, the middle of the 
C(N, M) band has to be used, cf. the dotted lines in the right column of Fig. 3. 

Using these interpolation curves to get the correction factors, the 1-D FM's can be 
corrected according to Eq.(8). The results are shown in Fig.4. The sawteeth seen in 
Fig. 2 are largely eliminated. 



4 Higher-dimensional Correction 

Having obtained the correction factors for one- dimensional FM's with non-integer M, 
let us turn to the correction of higher-dimensional FM's. 

First, consider the 2-D case. Let Mi, M 2 be the non-integer partitions in directions 
1 and 2, respectively, N ± , N 2 their integer parts. (M 1 — Ni + a±, M 2 = N 2 + a 2 , 
0<ai<l,0<a 2 <l). The definition of a non-integer FM, Eq.(8), is extended to 

1 Nl N2 (n (n - D) 
F 2 ( Nl , M i; N 2 , M 2 ) = - J- £ E { , U Z • ( U ) 

^ViiV 2 mi=lma=1 \n mi , m2 ) 

The average is taken over the area Ni ■ N 2 , shown as blank area in Fig.5a. Since the 
distribution of (n miifn2 (n miim2 — 1)) is in general not flat, a correction should be applied 
to take the shaded area in Fig. 5a into account. 

The one-dimensional correction factors Ri, R 2 of directions 1 and 2 can take care 
of the areas ASi = iV 2 (Mi - jVi) and AS 2 = N X (M 2 - N 2 ), respectively. However, the 
crossing area ASu = (Mi — iVi)(M 2 — N 2 ) has still to be considered. We evaluate it 
using an area counting method. Let the area ratios be 

, _ AS l2 _ M 1 -N 1 _ AS 12 _ M 2 -iV 2 

bl ~ a^-^vT" ; b2 ~ asT-^v, - • (12) 
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The correction factor due to ASu can be written as 



D _ p b 2/2 o 6 l/ 2 



Noting that i?i « 1, R 2 ~ 1, we have approximately 



i?i2 = l + -[&i( J R 2 -l) + 62(^1-1)]. 
The correction factor i?( 2 ) for 2-dimensional FM is then 



F 2 (M 1 ,M 2 ) 



— Ri ■ R 2 ■ -R12, 
1 



f 2 (n u m i; n 2: m 2 ). 



(13) 

(14) 
(15) 



In the left column of Fig. 6 are shown the 2-D FM's with real partition number M 
without correction. A strong "sawtooth" effect can be seen. Using correction factors 
obtained as above, Eq. (13-15), the large sawteeth disappear, cf. the right column of 
Fig.6. 

Similar considerations lead to the correction factor for the 3-dimensional case, cf. 
Fig.56: 



F 2 (M ± ,M 2 ,M 3 ) 



F 2 (N 1 ,M 1 ;N 2 ,M 2 ;N 3 ,M 3 ) 



R^ — Ri ■ R 2 ■ R 3 ■ R\ 2 ■ R 23 ■ i?3i • -R123, 



i?12 = l + -[&l(i?2-l) +62(^1-1)] 

R2z = l+ l -MRz-l) + h(R 2 -l)} 
R31 = l + ^[6 3 (i2i-l) + 6i(i23-l)] 



#123 = 1 + ^[bib 2 (R 3 - 1) + b 2 b 3 (R x - 1) + Mi (#2 - 1) 



61 = 
b 2 = 

h = 



(16) 
(17) 



AV 12 


AV 31 
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AV 3 
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The resulting 3D plots before and after correction are shown in Fig. 7. 
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5 Discussion 

In this paper we have developed a systematic method for an FM analysis with real 
(integer and non-integer) partition M. Correction factors are introduced in order to 
minimize the influence of the variation of f2(jn) = (n m (n m — 1))/ (n m ) 2 . The corrected 
results for non-integer M lie on smooth curves together with the integer M points. 

Note that only one-dimensional correction factors need to be extracted from the 
experimental data. In higher-dimensional correction the correction factors determined 
in one dimension have been used together with a simple geometrical consideration. This 
confirms the self-consistency of the method. 
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Figure captions 

Fig.l f2(m) = (n m (n m — 1))/ (n m ) 2 , the distribution in rapidity y transformed into the 
cumulant variable. 

Fig. 2 The one-dimensional plots of lni 7 ^ versus InM for real M (from M — 1 to 
M = 3.6) as defined in Eq.(5). The full circles correspond to integer M, while the open 
circles correspond to non-integer M. 

Fig. 3 The correction matrix C(N, M) as a function of N/M for M — 3 - 40. 

Fig. 4 The one-dimensional plots of lnF 2 versus InM for real M (from M — 1 to 

M = 3.6) as defined in Eq.(5) corrected according to Eq.(8). The full circles correspond 

to integer M, while the open circles correspond to non-integer M. 

Fig. 5 Schematic plot of the phase-space area (volumn) used in the average of a two- 

(three-) dimensional FM for non-integer partition M. The shaded areas (volumns) are 

those which have to be accounted for in the correction factors. 

Fig. 6 The two-dimensional plots of lni 7 ^ versus InM. The left column are the results 
of real M as defined in Eq.(5). The right column are the results after correction. 
Fig. 7 The three-dimensional plots of lnF 2 versus InM. The left column is the result of 
real M as defined in Eq.(5). The right column is the result after correction. 
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